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ABSTRACT 

We investigate numerically the ability of three models (jet, structured outflow, energy injection) to 
accommodate the optical light-curve breaks observed in 10 GRB afterglows (980519, 990123, 990510, 
991216, 000301c, 000926, 010222, 011211, 020813, and 030226), as well as the relative intensities 
of the radio, optical, and X-ray emissions of these afterglows. We find that the jet and structured 
outflow models fare much better than energy injection model in accommodating the multiwavelength 
data of the above 10 afterglows. For the first two models, a uniform circumburst medium provides a 
better fit to the optical light-curve break than a wind-like medium with a r~ 2 stratification. However, 
in the only two cases where the energy injection model may be at work, a wind medium is favoured 
(an energy injection is also possible in a third case, the afterglow 970508, whose optical emission 
exhibited a sharp rise but not a steepening decay). The best fit parameters obtained with the jet 
model indicate an outflow energy of 2 — 6 X 10 50 ergs and a jet opening of 2° — 3°. Structured outflows 
with a quasi-uniform core have a core angular size of 0.7° — 1.0° and an energy per solid angle of 
0.5 — 3 X 10 53 ergs sr — 1 , surrounded by an envelope where this energy falls-off roughly as 9~ 2 with 
angle from the outflow axis, requiring thus the same energy budget as jets. Circumburst densities 
are found to be typically in the range 0.1 — 1cm -3 , for either model. We also find that the reverse 
shock emission resulting from the injection of ejecta into the decelerating blast wave at about 1 day 
after the burst can explain the slowly decaying radio light-curves observed for the afterglows 990123, 
991216, and 010222. 

Key words: gamma-rays: bursts - ISM: jets and outflows - radiation mechanisms: 
non-thermal - shock waves 



1 INTRODUCTION 

Since the prediction of radio (Paczyhski & Rhoads 1993) and 
optical transients (Meszaros & Rees 1997) associated with 
GRBs, more than 100 afterglows have been observed (in- 
cluding X-ray transients). Good monitoring in all three fre- 
quency domains has been achieved for about 20 afterglows; 
the radio, optical, and X-ray flux was observed to decay as a 
power- law, F v oc t~ a (a > 0), confirming the expectations. 
For 10 of these well-observed afterglows - 980519, 990123, 
990510, 991216, 000301c, 000926, 010222, 012111, 020813, 
030226 - the optical light-curve decay exhibits a steepening 
at about 1 day after the burst, with an increase Aa in the 
temporal index as low as 0.4 (afterglow 991216) and as high 
as 1.8 (afterglow 000301c). Optical light-curve breaks have 
or may have been observed in other afterglows; they are not 
included in the sample of afterglows modelled in this work 
because those afterglows have been adequately monitored at 
only one optical frequency. 

With the exception of the afterglow 010222, for which 
one, late (10 days) measurement suggests that a break oc- 
curred also in the X-ray light-curve, currently available X- 



ray observations do not extend sufficiently before and after 
the optical break time to prove that the light-curve break is 
achromatic. Within the measurement uncertainties, a single 
power-law fits well the decay of all adequately monitored X- 
ray afterglows: 990123, 990510, 991216, 000926 (and 010222 
until the last measurement). 

Breaks, in the form of peaks, have been observed in the 
radio emission of all 10 afterglows above, usually at ~ 10 
days after the burst, i.e. often occurring after the optical 
light-curve break. Given that the radio and optical breaks 
are not simultaneous, they cannot both arise from the dy- 
namics of the afterglow, the structure of the GRB ejecta, or 
some property of the circumburst medium (CBM), as these 
mechanisms should yield achromatic breaks. Furthermore, 
given that, before their respective breaks, the radio emis- 
sion rises while the optical falls-off, these temporal features 
cannot be due to the passage of the same afterglow contin- 
uum break through the observing band. 

With the possible exception of the afterglow 990123, 
there is no evidence for an evolution (softening) of the op- 
tical continuum across the light-curve break for the above 
10 afterglows (e.g. Panaitescu 2005). Furthermore, the only 
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spectral feature whose passage could yield the large steep- 
ening Aa observed in some cases - the peak of the forward 
shock (FS) continuum*- should yields a rising or flat optical 
light-curve before the break, contrary to what is observed. 
Therefore, it is the radio peak which should be attributed to 
a spectral break crossing the observing domain. That spec- 
tral break should be the injection frequency Vi, which, for 
reasonable shock microphysical parameters, should reach the 
radio at around 10 days. Indeed, for the afterglow 991208, 
there is observational evidence (Galama et al. 2000) that 
the peak of the afterglow continuum crosses the 10-100 GHz 
domain at that time. Further evidence for the passage of a 
spectral break through radio is provided by that the peak 
time of the radio flux of the afterglow 030329 increases with 
decreasing observing frequency (Frail et al. 2005). 

Today, the generally accepted reason for the optical 
light-curve break is the narrow collimation of the GRB 
ejecta. As predicted by Rhoads (1999), if the GRB ejecta 
are collimated, then the afterglow light-curve should exhibit 
a steepening when the jet begins to expand sideways. More 
than half of the steepening Aa is due by the finite angu- 
lar opening of the jet: as the GRB remnant is decelerated 
(by sweeping-up of the CBM) and the relativistic Doppler 
beaming of the afterglow emission decreases, an ever increas- 
ing fraction of the emitting surface becomes visible to the 
observer; when the jet edge is seen, that fraction cannot 
increase any longer and the afterglow emission exhibits a 
faster decay (Panaitescu, Meszaros & Rees 1998). Because 
it arises from the blast wave dynamics, a light-curve break 
should also be present at radio wavelengths at the same time 
as the optical break. However, radio observations before 1 
day are very scarce and strongly affected by Galactic inter- 
stellar scintillation (Goodman 1997) until after 10 days, thus 
they cannot disentangle the jet-break from that arising from 
the passage of the FS peak frequency. 

Another mechanism for the optical light-curve breaks 
has been proposed by Rossi, Lazzati & Rees (2002) and 
Zhang & Meszaros (2002): if GRB outflows are endowed 
with an angular structure (i.e. non-uniform distribution of 
the ejecta kinetic energy with direction), as first proposed 
by Meszaros, Rees & Wijers (1998), then a steepening of the 
afterglow decay would arise when the brighter, outflow sym- 
metry axis becomes visible to the observer. In this model, 
the stronger the angular structure is, the larger the break 
magnitude Aa should be. 

A third mechanism for breaks rests on the proposal of 
Paczyhski (1998) and Rees & Meszaros (1998) that the FS 
energizing the CBM could be refreshed by the injection of a 
substantial energy through some delayed ejecta which were 
released at the same time with the GRB-producing ejecta, 
but had a smaller Lorentz factor, or were ejected sometime 
later, and which catch up with the decelerating FS during 
the afterglow phase. Fox et al. (2003) have proposed that the 
early (0.003-0.1 day), slow decay of the optical emission of 
the afterglow 021004 is caused by such an injection process. 
In this scenario, when the energy injection episode ends, the 



* This is the smallest of the synchrotron frequency Vi correspond- 
ing to the typical post-shock electron energy (which we call "in- 
jection frequency") and to the electrons which cool radiatively on 
a dynamical timcscale (the cooling frequency u c ) 



FS deceleration becomes faster and the afterglow emission 
should exhibit a steepening. 

Note the various origins of the afterglow light-curve 
break in each model. In the Jet model, the break is caused 
by the changing outflow dynamics when the jet starts to 
spread and by the outflow's geometry. In the Structured 
Outflow (SO) model the origin is, evidently, the outflow's 
anisotropic surface brightness. For both these models, spe- 
cial relativity effects play an important part. In the Energy 
Injection (EI) model, the break originates in the altered 
outflow dynamics at the time when the energy injection sub- 
sides. 

The purpose of this work is to compare the ability of 
these three models in accommodating 

i) the shape of the light-curve breaks observed in the optical 
emission of the afterglows 980519, 990123, 990510, 991216, 
000301c, 000926, 010222, 012111, 020813, 030226 and 

ii) the relative intensity of the radio, optical, and X-ray 
emissions of these afterglows, 

for either a uniform (i.e. homogeneous) CBM, or one with 
a r~ 2 density radial stratification, as expected if GRBs pro- 
genitors are massive stars (Woosley 1993, Paczyhski 1998). 

For the first task above, we performed an analytical 
test of the models (see table 2 of Panaitescu 2005), based 
on comparing 

i) the pre- and post-break optical light-curve decay indices, 
ai and CV2, 

ii) those in the X-rays, 

Hi) the slopes (5 and (3 X of the power-law optical and X-ray 
continua {F a ocf^ 13 " and F x oc v~ !)x ), and 
iv) the optical-to-X-ray spectral energy distribution (SED) 
slope, (3ox = ln(F x /F )/ln(v x /v ), 

with the relations among them expected for each model. 
We note that, in the framework of the relativistic fireballs 
(Meszaros & Rees 1997), the optical and X-ray decay indices 
and SED slopes are tightly connected, as only one continuum 
feature, the cooling frequency v c (Sari, Narayan & Piran 
1998), can be between these domains at the times when 
observations were usually made (0.1-100 days), and that the 
SED slope increases by a fixed amount, <5/3 = 1/2, across this 
spectral break. 

Including the radio afterglow emission in the analyti- 
cal test of the three break models is less feasible and often 
unconstraining. First, that the afterglow radio flux is mod- 
ulated by diffractive and refractive interstellar scintillation 
makes it difficult to determine accurately the radio SED 
slope (see figs. 4 and 5 of Frail, Waxman & Kulkarni 2000a 
for the best monitored radio afterglow - 970508). Second, 
after the injection frequency has fallen below the radio (i.e. 
during the decay phase of the radio light-curve), we do not 
expect, in general, any spectral break to be between radio 
and optical, hence the radio and optical light-curve indices 
should be the same. This is, indeed, the case for most af- 
terglows; nevertheless, there are a few troubling exceptions: 
over 1-2 decades in time, the radio emission of the after- 
glows 991208, 991216, 000926, and 010222 exhibits a much 
shallower decay than that observed at optical wavelengths 
at the same time or prior to the radio decay. An analytical 
investigation of the various possible ways to decouple the 
radio and optical light-curves (Panaitescu & Kumar 2004) 
has led to the conclusion that, the anomalous radio decay 
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is due to a contribution from the reverse shock to the radio 
emission. 

However, radio observations provide an indirect con- 
straint on them because the flux and epoch of the radio peak 
determine the FS synchrotron peak flux, injection frequency, 
and self-absorption frequency. These three spectral proper- 
ties constrain the afterglow parameters pertaining to the 
outflow dynamics (energy per solid angle, jet initial open- 
ing, medium density) and emission (magnetic field strength, 
typical post-shock electron energy), i.e. parameters which 
determine the shape and epoch of the optical light-curve 
break, as well as the location of the cooling frequency. The 
best way to take into account the constraints arising from 
the radio emission and to test fully the three break mod- 
els is to calculate numerically the afterglow dynamics and 
emission, and to fit all the available measurements. Data 
fitting also allows the determination of the various model 
parameters. 



2 MODELS DESCRIPTION 

The basic equations employed in our numerical calcula- 
tions of the outflow dynamics are presented in Panaitescu & 
Kumar (2000) for spherical outflows, Panaitescu & Kumar 
(2001) for the Jet model, Panaitescu & Kumar (2003) for the 
SO model (structured outflows), and Panaitescu, Meszaros 
& Rees (1998) for the EI model (energy injection). The ba- 
sic equations for the calculation of the afterglow spectral 
features (synchrotron and inverse Compton peak fluxes, ab- 
sorption, injection, and cooling frequencies), and of the af- 
terglow emission at any wavelength are given in Panaitescu 
& Kumar (2000, 2001). Equations for the spectral character- 
istics can be also found in Sari et al. (1998), for the dynamics 
and emission from spherical blast-waves interacting with a 
uniform medium in Waxman, Kulkarni & Frail (1998), Gra- 
not, Piran & Sari (1999), Wijers & Galama (1999) and for 
a wind medium in Chevalier & Li (2000). The dynamics of 
blast-waves with energy injection is also treated in Sari & 
Meszaros (2000). Rhoads (1999) provided a detailed treat- 
ment of the dynamics and emission from jets interacting 
with uniform media (see also Sari, Piran & Halpern 1999). 



assumed to increase at the co-moving sound speed, and the 
kinetic energy per solid angle is approximated as uniform 
during the spreading. Radiative (synchrotron and inverse 
Compton) losses are calculated from the electron distribu- 
tion and magnetic field strength. 

In the Jet model, the light-curve break is seen when the 
jet edge becomes visible to the observer. For an observer 
located close to the jet axis, this time can be approximated 
by T(tj e t) = djet, using the initial jet opening and ignoring 
the lateral spreading that occurred until the jet-break time. 
If the FS dynamics were adiabatic, then energy conservation 
would lead to 

r(t) = 8.6 0Eo,5 3 /no) 1/8 [W(l + ^)]r 3/8 (2) 
for a uniform medium and 

r(t) = 15 (£ ,53M*) 1/4 [td/(l + z)])- 1/A (3) 

for a wind, where td is the observer time in days. Then the 
jet-break time is given by 

t jet ~ 0.7{z + l)(£o,5 3 no ^fet.-i) 173 d (4) 
for a uniform CBM and 

t jet ^ 5( Z + l)£ ,53^;:X 4 e*,-l d (5) 

for a wind. Given that observer locations off the jet axis 
(but within 9 jet) have a small effect on the resulting after- 
glow light-curves (Granot et al. 2002), we consider that, in 
the Jet model, the observer is always on the jet axis. Fur- 
thermore, we ignore the possible existence of a counter-jet 
whose emission would become visible to the observer when 
the semi-relativistic dynamics sets, i.e. it would affect the 
radio afterglow emission beyond 100 days after the burst. 

At the jet-break time, the afterglow decay index at ob- 
serving frequency v steepens from that corresponding to a 
spherical outflow (Sari et al. 1998): 

1 r 3p - 3 , v < v c h unif CBM 
Qi = - • < 3p - 2 , v c < v & any CBM , (6) 
I 3p - 1 , v < v c & wind CBM 

to Q2 = p (Rhoads 1999). In equation (6), p > is the power- 
law index of the post-shock electron distribution with energy 
(equation [16]). 



2.1 Jets 

The jet dynamics is determined by the initial jet opening 
9 j e t, the ejecta initial kinetic energy per solid angle Eq (or, 
equivalently, by the jet energy Ej e t = n9 2 j et Eo), and the 
CBM particle density n or, in the case of a wind of a massive 
star, by the ratio of the mass-loss rate to the wind speed, 
(dM/dt)/v. For convenience, the latter parameter is given 
normalized to a mass-loss rate of 10~ 5 M yr _1 at a speed 
of 10 3 kms _1 , i.e. for a wind 

n(r) = 0.3 A»r w 2 cm -3 (1) 

where r is the FS radius, and we used the usual notation 
Q„ = Q(cgs units)/ 10™. 

The deceleration of the jet, as it sweeps-up the CBM, 
is calculated assuming that the post-shock gas has the same 
internal energy per mass (i.e. temperature) as that imme- 
diately behind the FS, which is equal to the bulk Lorentz 
factor of the shocked CBM, F. The lateral size of the jet is 



2.2 Structured Outflows 

The dynamics of structured outflows is calculated similarly 
to that of jets. To track the lateral fluid flow and the change 
of the kinetic energy per solid angle, the outflow surface is 
divided in infinitesimal rings and we consider that each ring 
spreads at a rate proportional to the local sound speed and 
to the ring width. This prescription for lateral flow reduces 
to that given above for a jet if the outflow is uniform and 
collimated. 

In this work we consider only axially symmetric power- 
law outflows, whose angular distribution of the ejecta kinetic 
energy per solid angle E is given by 

E{9 < 9 core ) = E , E(9 > 9 core ) = E (9/9 core )- q , (7) 

with q > 0. The uniform core is used to avoid a diverging 
outflow energy for q > 2. This power-law structure is suffi- 
ciently complex, given the constraining power of the avail- 
able data. Given that, in this model, the light-curve break 
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is due to the brighter core becoming visible, the observer 
location 9 bs relative to the symmetry axis is a crucial pa- 
rameter in determining the break time (Rossi et al. 2002), 
which is the same as given in equations [4] and [5] but with 
B obs in place of 9 jet - 

Analytical results for the light-curve pre- and post- 
break decay indices can be obtained if the observer is located 
within the core (9 bs < 9 core), but only numerically for outer 
locations^. In the former case, the light-curve break occurs 
when the edge of the uniform core becomes visible to the ob- 
server, and the light-curve decay index at a frequency v > Vi 
increases from that given in equation (6) to 



Q.2 



3p- 
3p- 



3 + |g v < v c 
2 + q v c < v 



for a homogeneous medium and 

_J_ f3p-l- lg(p-l) v<V c 

2 4-q \3p-2- ig(p-2) v c < v ' 



(8) 



(9) 



for a wind medium (Panaitescu & Kumar 2003), provided 
that the structural parameter q < q, where, for a uniform 
medium, q — 8/(p + 4) if v < v c and q — 8/(p + 3) if v > v c . 
For a wind, the two values of q are swapped. For q < q, the 
emission from the outflow envelope (9 > 9 core ) is brighter 
than that from the core, and sets the light-curve post-break 
decay index. For q > q, the core dominates the afterglow 
emission and the post-break decay index is 



Ct2 



ai + 



3/4 
1/2 



unif CBM 
wind CBM 



(10) 



i.e. Aa is just the steepening produced by the finite angular 
extent of the emitting outflow. In this case the structure 
outflow model reduces to a that of a jet with sharp edges 
and no sideways expansion. 

For observers located outside the core, the break mag- 
nitude Aa can be larger (see figs. 2, 3, and 4 in Panaitescu 
& Kumar 2003). 



2.3 Energy Injection 

The injection of energy in the FS has two effects: it miti- 
gates the FS deceleration and generates a reverse shock (RS) 
which energizes the incoming ejecta and contributes to the 
radio afterglow emission. 

For ease of interpretation, we consider an energy injec- 
tion which is a power-law in the observer time: 



dEi/dt{t < t ff) oc t e 1 , dE,/dt{t > t ff) = , 



(11) 



where dEi/dt is the rate of the influx of energy per solid 
angle. Evidently, this injection has an effect on the afterglow 
dynamics only if the total added energy, Ei, is comparable 
or larger with that existing in the afterglow after the GRB 
phase, Eo- In this case, the light-curve decay index during 
the injection process is 



t For this reason, in our previous analytical assessment of the 
three models (Panaitescu 2005), we have restricted our attention 
only to the o bs < Score case, which yields breaks with a smaller 
magnitude Aa 



ai = i(3-e)p-i(l+e)- 



3 , v < v c & unif CBM 

2 , v c < v & any CBM ,(12) 

1 , v < v c & wind CBM 



(Panaitescu 2005). After the energy injection subsides, the 
post-break index a? has the value given in equation (6) (for 
an adiabatic, spherical blast- wave). 

Numerically, the energy injection is modelled by first 
calculating the dissipated energy and that added as kinetic 
energy during the collision between an infinitesimal shell 
of delayed ejecta and the FS. The partition of the incom- 
ing ejecta energy is determined from energy and momentum 
conservation, and depends only on the ratio Ti/r, with Fi 
being the Lorentz factor of the incoming ejecta. Numerically, 
this factor is obtained from the kinematics of the catching-up 
between the freely flowing delayed ejecta and the decelerat- 
ing FS. It can also be calculated analytically, if it is assumed 
that all the ejecta were released on a timescale much shorter 
than the observer time when the catching up occurs. With 
this assumption, it can be shown that the ejecta-FS contrast 
Lorentz factor is constant: 



Fi/F = (1 + e) 



-1/2 



2 

V2 



unif CBM 
wind CBM 



(13) 



Once the dissipated fraction is known, we track the adiabatic 
conversion of the internal energy into kinetic, as described in 
Panaitescu et al. (1998). Although the concept of adiabatic 
losses implies that there will be a dispersion in the Lorentz 
factor of the swept-up CBM, we ignore it and, just as for 
the other two models, we assume that all the fluid behind 
the FS moves at the same Lorentz factor T. 

Even if the injected energy is negligible, the RS crossing 
the delayed ejecta can be of relevance for the lower frequency 
afterglow emission at days after the burst. To calculate this 
emission, we set the typical electron energy behind the RS 
by assuming that all the dissipated energy is in the shocked 
ejecta, which would be strictly correct only if the density 
of the shocked CBM were much larger than that of the de- 
layed ejecta. Otherwise, this assumption leads to an over- 
estimation of the post-RS electron energy. After the energy 
injection ceases, we track the evolution of the RS electron 
distribution, subject to adiabatic and radiative cooling. The 
tracked electron distribution is used for the calculation of the 
RS continuum break- frequencies. In contrast, for all three 
models, the FS electron distribution is set using the current 
Lorentz factor T of the FS (see below). This last approxima- 
tion is more appropriate for a uniform medium than for a 
wind, as in the former case dM/dr oc r 2 while for the latter 
dM/dr = constant. 



2.4 Model Parameters 

As presented above, the Jet model has two dynamical pa- 
rameters: the initial kinetic energy per solid angle Eo and 
the initial jet opening 9 jet- The SO model has three such pa- 
rameters: the energy per solid angle Eo on the axis (and in 
the uniform core), the core angular size 9 core, and the struc- 
tural parameter q for the outflow envelope. The EI model 
has three dynamical parameters: the total injected energy 
Ei, the temporal index e of the injection law, and the time 
t ff when the injection episode ends, the initial Eo being less 
relevant if Eo <^ Ei . In addition, all models have another 
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parameter which determines the outflow dynamics: the ex- 
ternal medium density n (or the wind parameter A*). The 
observer location 6 obs relative to the outflow symmetry axis 
is relevant only for the SO model, but not for the Jet model, 
as long as o bs < 6 jet , or for the EI model, where the outflow 
is spherically symmetric. 

The calculation of the synchrotron and inverse Comp- 
ton emissions (and of the radiative losses), requires three 
more parameters: one for the magnetic field strength, pa- 
rameterized by the fraction eb of the post-shock energy in 
it: 

B 2 /(Sty) = n (4r + 3)(r - 1) m p c ■ eb , (14) 

n being the co-moving frame particle density in the un- 
shocked fluid, m p the proton mass, one for the minimal elec- 
tron energy behind the shock, parameterized by the fraction 
Ei of the post-shock energy in electrons if all had the same 
energy 7, m e c 2 : 

7; m e = (r — l)m p ■ Ei , (15) 

7i being the electron random Lorentz factor and m e the 
electron mass, and the index p of the electron power-law 
distribution with energy: 

dN/de oc e~ v . (16) 

Equations (14) and (15) apply also to the RS if Y is replaced 
by the Lorentz factor of the shocked ejecta as measured in 
the frame of the incoming (unshocked) ejecta. The parame- 
ters Ei and eb determine the synchrotron characteristic fre- 
quency Vi corresponding to the typical electron energy, as 
well as the self-absorption and cooling frequencies. 

We assume that the three microphysical parameters e b , 
Ei, and p, have the same value behind both the forward and 
reverse shock. 

Summarizing, the Jet model has six free parameters, 
the EI model seven, and the SO model eight. The V-band 
extinction due to dust in the host galaxy, Ay, which affects 
the observed slope of the optical SED and the overall optical 
flux, is an extra parameter for all models. To determine the 
host frame Ay extinction, we assume an SMC-like reddening 
curve $ for the host galaxy and use the parameters inferred 
by Pei (1992) for SMC. 

All these parameters are constrained by the multiwave- 
length afterglow data in the following way. The pre-break 
decay index a\ of the optical light-curve determines the elec- 
tron distribution index p for the Jet and SO models (eqs.[6], 
[8], [9]), and the energy injection index e for the EI model 
(eq.[12]). The post-break decay index 02 overconstrains the 
electron index p for the Jet model, determines the struc- 
tural parameter q for the SO model, and the electron in- 
dex p for the EI model§. The jet opening 9j e t for the Jet 
model and the observer location 9 obs for the SO model (if the 

t In principle, if a simple power-law A„ oc u K reddening curve is 
assumed, then numerical fits to the afterglow optical continuum 
and the X-ray flux could constrain both Ay and the index n 
§ Therefore, for the SO and EI models, the double constraint on 
the electron index p specific to the Jet model is relaxed. However 
this does not guarantee that the fits obtained with the SO and 
EI models will always be better than those resulting for the Jet 
model, as there are differences in the blast- wave dynamics of these 



observer is outside the core) or the core angular size 6 CO re 
(if the observer is within the core) are set by the epoch of 
the optical light-curve break (eqs.[4] and [5]). The same ob- 
servable determines the time t a f / when the energy injection 
ceases in the EI model. For the SO model, whenever the 
observer is outside the outflow core, the 9 core is not well- 
constrained because it has a weak effect on the resulting 
afterglow emission. The outflow energy Eo (or the injected 
one Ei), the medium density n (or the wind parameter A-,), 
and the two microphysical parameters e< and eb are con- 
strained by the afterglow spectral parameters - flux at the 
peak of the spectrum, self-absorption, injection, and cooling 
frequencies - which in turn are constrained by the relative 
intensities of the radio, optical and X-ray emissions (e.g. 
Wijers & Galama 1999, Granot et al. 1999). That the self- 
absorption frequency is not always constrained by the radio 
observations accounts for part of the resulting parameter un- 
certainties, but note that the parameters Eo, n, eb, and Ei 
are also constrained by matching the radio emission at late 
times, when the blast-wave is only mildly relativistic, i.e. 
after the self-similar relativistic dynamics phase employed 
in analytical calculations of afterglow parameters. For the 
EI model, the initial outflow energy Eo is irrelevant, as the 
blast-wave kinetic energy is that injected starting at a time 
before the first observation. Finally, the host extinction Ay 
is constrained by the ratio of the observed optical flux to 
the intrinsic optical flux which is consistent with the X-ray 
emission (for the other afterglow parameters) . 



3 RESULTS OF NUMERICAL MODELLING 

Table 1 lists the reduced chi-square, xti f° r the Des t fits ob- 
tained with the Jet, SO, and EI models, for the two types of 
circumburst media: homogeneous and wind-like. For a uni- 
form CBM, the best fits with the EI model are very poor 
and, for almost all afterglows, they are substantially worse 
than those for a wind. This is so because, for a uniform 
CBM, the observed post-break optical decay index «2 (given 
by equation [6]) requires a higher electron index p than for a 
wind, leading to an optical SED slope (/?<,) and an optical- 
to-X-ray spectral index (f3 ox ) larger than those observed. 
Also, the EI model has difficulties in accommodating both 
the radio and X-ray fluxes, with the FS emission either over- 
estimating the observed radio flux or underestimating the 
X-ray flux. 

Tables 2, 3, 4, and 5 list the parameters of the best 
fits with xl < 4 for the Jet, SO, EI, and Jet+EI models, 
respectively, some of which are shown in Figs. 1-17. From 
the parameter uncertainties determined by us (Panaitescu 
& Kumar 2002) for the Jet model from the x 2_va riation 
around its minimum for eight GRB afterglows, we estimate 
the following uncertainties for the best fit parameters given 
in Tables 2-5: <r(lg£b) = ff(lg-Ei) = 0.3, ff(lgn) = 0.7, 
ff(lg A,) = 0.3, a(6 jet ) = 0.30,-et, cr((9 core ) = 0.2°, cr(lge s ) = 
1, o(£i) = 0.3e<, a(p) = 0.1, a(q) = 0.3, cr(e) = 0.2, a(A v ) = 
0.2 Ay. 

Figs. 1-17 display the best fits obtained with the 
Jet model for the afterglows 010222, 011211, 020813, and 

three models and specific limits on the post-break light-curve de- 
cay indices which they yield. 
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Table 1. Reduced chi-square xi of the best fits obtained for 10 GRB afterglows with three models (SO = structured outflow, EI = 
energy injection) for light-curve breaks and for two type of circumburst media: uniform (n = const) and wind-like (n oc r~ 2 ) 



/~< r_> 15 




l\ 


Jet 


Jet 


bU 


bU 






(1) 


(2) 


n = const 


n oc r~ 2 


n = const 


n oc r~ 2 


n oc r~ 2 


980519 


0.56 ±0.37 


73 


2.6 


1.8 


2.4 


1.4 


3.0 


990123 


0.48 ±0.21 


112 


2.0 


2.3 


1.8 


2.2 


1.5 


990510 


1.29 ±0.10 


101 


0.78 


3.1 


2.1 


4.6 


3.0 


991216 


0.40 ±0.20 


84 


2.0 


1.8 


1.2 


1.2 


3.4 


000301c 


1.83 ±0.18 


111 


4.4 


8.3 


3.3 


7.1 


10 


000926 


0.64 ±0.13 


145 


2.2 


3.5 


2.2 


2.8 


3.3 


010222 


0.88 ±0.08 


175 


2.2 


3.9 


1.7 


4.0 


4.7 


011211 


1.70 ±0.25 


88 


4.7 


8.7 


2.3 


4.7 


7.6 


020813 


0.54 ±0.06 


105 


1.6 


2.6 


1.1 


2.6 


1.6 


030226 


1.48 ±0.07 


112 


8.5 


17 


4.0 


11 


16 



(1) : magnitude of optical light-curve break, defined as the increase of the exponent of the power-law flux decay 

(2) : number of radio, optical, and X-ray measurements 

Table 2. Jet model parameters for the best fits of Table 1 with \i < 4 
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0.1 




0.8 
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fig.2 


in PK02 






2.0 


0.2 


6.7 


-0.7 


-1.4 


2.4 


0.06 






990123 


0.8 




0.9 


2.2 


-3.8 


-2.7 


1.5 


0.18 




4 






0.2 


0.7 


2.3 


-2.8 


-3.1 


1.5 


0.16 






990510 


0.3 




0.2 


3.1 


-2.3 


-1.6 


1.8 


0.04 


fig.3 


in PK01 






0.4 


1.5 


2.1 


-3.2 


-1.6 


1.8 


0.05 






991216 


0.04 




0.6 


1.5 


-1.5 


-2.5 


1.7 


0.05 




6 






0.2 


0.8 


2.6 


-2.3 


-1.9 


1.7 


0.05 






000926 


22 




0.1 
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-1.3 


-1.1 


2.4 


0.14 


fig.5 


in PK02 
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0.2 


5.1 


-1.3 


-1.1 


2.4 


0.12 






010222 


0.09 




0.8 


2.5 


-3.1 


-2.0 


1.8 


0.28 




10 






0.1 


0.2 


3.8 


-2.1 


-2.5 


1.5 


0.30 






020813 


0.07 




0.7 


2.3 


-2.9 


-2.1 


1.9 


0.25 




14 






0.2 


0.4 


4.8 


-3.5 


-2.3 


1.5 


0.29 







Table 3. Structured outflow parameters for the best fits of Table 1 with \1 < 4 
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( cm -3 ) 
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8 obs 1 8 core 
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Ay 


Fig 


980519 


4 




0.6 


0.9 


1.4 


0.9 


-3.4 


-1.5 


2.7 


0.07 


1 






0.6 


1.0 


1.0 


2.3 


1.6 


-3.8 


-1.4 


2.6 


0.07 




990123 


1 




0.5 


0.7 


1.7 


0.7 


-3.2 


-2.2 


1.6 


0.15 








0.1 


0.2 


1.1 


1.5 


1.3 


-3.0 


-2.6 


1.5 


0.14 




990510 


2 




0.7 


0.7 


1.8 


3.0 


-2.7 


-1.4 


1.8 







991216 


0.7 




3 


0.5 


1.6 


1.5 


-3.9 


-2.0 


1.6 











0.4 


3 


0.5 


1.2 


0.5 


-3.8 


-1.7 


1.7 


0.05 




000301c 


0.2 




2.5 


0.8 


2.3 


3.3 


-2.6 


-1.6 


2.4 





8 


000926 


3 
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0.7 


2.6 


3.4 


-2.8 


-1.3 


2.4 


0.14 
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0.5 
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0.9 


2.6 


1.9 


-2.7 


-1.5 


2.6 


0.18 




010222 


0.09 




0.6 


1.7 


2.1 


0.7 


-3.7 


-1.9 


1.8 


0.27 
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1.6 
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0.11 
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0.9 
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1.1 


2.2 
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-3.3 


-1.3 


2.3 
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0.06 
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0.20 








0.06 
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1.0 


2.0 


-2.6 
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0.20 
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1 




2 


0.5 


2.8 


3.6 


-3.6 


-1.5 


2.4 


0.03 


17 



© 2005 RAS, MNRAS 000, 000-000 



Jets, Structured Outflows, and Energy Injection in GRB Afterglows 

Table 4. Best fit parameters obtained with the Energy Injection model 
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990123 


0.003 


0.08 


0.7 
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-2.1 


-1.0 


2.3 


0.04 


5 


020813 


0.13 


0.7 


0.3 
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Table 5. Best fit parameters obtained with the Jet+EI model 
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6.0 
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10- 3 
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0.8 


3.6 


0.2 
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-1.5 


2.0 







010222 


1.6 


10- 4 




0.7 


2.0 


0.06 


2.5 


-1.8 


-1.3 


2.3 


0.19 


11 




3.6 




0.06 


0.4 


5.0 


0.6 


1.7 


-3.0 


-1.7 


1.7 


0.19 





030226, which we did not present previously, and compara- 
ble or better fits obtained with the SO and EI models. With 
the exception of the Jet model and a uniform CBM for the 
afterglow 990510, all best fits presented are not acceptable in 
a statistical sense, as they have xl > 1- I n general, fits with 
xt < 4 appear adequate upon visual inspection, a significant 
fraction of the \ 2 arising often from optical measurements 
(which have the smallest uncertainties), indicating either a 
small, short-lived afterglow fluctuation or, perhaps, an oc- 
casional outlier resulting from underestimated observational 
errors. In some cases there are also systematic discrepancies 
between the model fluxes and observations, most often in the 
radio. As a general rule, if the fits are poorer than \ 2 > 4 
for both types of circumburst media then the fits are not 
shown, with the exception of those obtained with the Jet 
and SO models for the afterglows 011211 and 030226, for 
which we want to show that the latter model is much better 
than the former in accommodating the 1 day, sharp optical 
light-curve break. Also as a general rule, if the fits obtained 
with the Jet and SO models are comparable, only one is 
shown in the figures. 

For clarity, the figures contain only the radio optical 
frequencies where observations span the longest time, but 
the fitted data set included few or several other bands in 
these two domains (see figure captions). With the exception 
of the afterglow 000926, the available X-ray measurements 
are in only one band, which we use to infer the X-ray flux at 
the mid energy of that band. The amplitude of interstellar 
scintillation, whose calculation is based on the treatment 
and maps given by Walker (1998), is indicated with vertical 
bars at the time of the radio observations. In all figures, the 
best fit for a uniform CBM is shown with continuous lines 
while that for a wind medium with dotted lines. 



3.1 980519 

The best fit obtained with the Jet model is discussed in 
Panaitescu & Kumar (2002). The addition of more radio 
data and some optical outliers (which were previously ex- 
cluded) lead now to a worse fit (x -wise) but to the same 



afterglow parameters, within their uncertainty. The best fit 
for the SO model, shown in Fig. 1, has a xl comparable to 
that for the Jet model (Table 1) but is qualitatively poorer, 
as it overestimates the observed early radio. Although the 
measured radio fluxes are within the amplitude of the fluc- 
tuations caused by the interstellar Galactic gas, such a con- 
stant offset cannot be caused by scintillation (unfortunately, 
X 2 -statistics does not reflect the smaller probability of sys- 
tematic discrepancies) . 



3.2 990123 

The best fit obtained with the Jet model and a uniform 
medium is presented in Panaitescu & Kumar (2001). There 
we assumed that the radio flare seen at 1 day (Kulkarni et 
al. 1999) arises in the GRB ejecta energized by the RS, and 
included in the modelling only the radio measurements at 3- 
30 days, when the 8 GHz flux of this afterglow is less than 
40/iJy. This sets an upper limit on the FS peak flux, F p , 
which, together with that F p oc n 1 / 2 before the jet-break 
and F p oc n 1 / 6 after that, requires a small CBM density: 
n < I0- 2 cm -3 . 

Including the radio measurements before a few days, 
as well as the millimeter data, the K-b&nd fluxes after 10 
days (with which the if-band light-curve appears to be a 
single power-law, unlike the i?-band emission, which exhibits 
a steepening at ~ 2 days), and other optical measurements 
previously left out of our modelling, we find now a poorer 
fit (xl = 2.4 for a uniform CBM, xl = 5.1 for a wind), but 
with the same afterglow parameters as before. 

We have also tested the ability of the Jet model includ- 
ing the emission from the RS crossing the GRB ejecta to 
explain the optical flash observed at 20-600 seconds (Ak- 
erlof et al. 1999), which peaked 50 s after the burst at about 
0.8 Jy, and the 1 day radio flare, as proposed by Sari & 
Piran (1999). Given that the observer injection timescale 
(~ 100 s) is comparable to the burst duration, the approx- 
imation of instantaneous ejecta release and its consequence 
given in equation (13) may be invalid. For this reason, we let 
the ejecta Lorentz factor, T;, to be a free parameter of the 
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fit. Together with the injected energy, Ei, it determines the 
number of the radiating electrons in the GRB ejecta and the 
Lorentz factor of the RS, i.e. the peak flux and frequency of 
the synchrotron RS emission. 

The best fits obtained with the Jet model interacting 
with a tenuous medium (so that the FS emission accommo- 
dates the late radio measurements) and emission from the 
GRB ejecta are shown in Fig. 2. We find that, for the same 
microphysical parameters s, and eb behind both shocks, the 
peak RS optical emission is dimmer than observed by a fac- 
tor 50 for a uniform CBM and a factor 5 for a wind, the 
latter yielding a good fit to the remainder of the early op- 
tical measurements (but not as good as a uniform medium 
to the 1 day optical break). We also find that, under the 
same assumption of equal microphysical parameters, the RS 
emission peaks at 8 GHz at 0.1 day, i.e. a factor 10 too 
early. In other words, taking into account the cooling of the 
RS electrons accelerated at ~ 100 s and the decrease on the 
synchrotron self-absorption optical thickness of the RS at 
radio frequencies, we do not find a set of afterglow parame- 
ters for which the RS emission can peak in the radio as late 
a 1 day while, for the same microphysical parameters, the 
FS emission can accommodate the rest of afterglow obser- 
vations. To explain the optical flash of the afterglow 990123, 
the GRB ejecta must have a larger parameter eb than for 
the forward shock, indicating that the ejecta is magnetized 
(Zhang, Kobayashi & Meszaros 2003, Panaitescu & Kumar 
2004), a feature that may be required for by the early op- 
tical emission of the afterglow 021211 as well (Kumar & 
Panaitescu 2003, Zhang et al. 2003). 

However, for the low density Jet model, an injection of 
ejecta in the FS at about the same time when the radio 
flare is seen allows this model to accommodate that flare, 
as shown in Fig. 3. Note that the injected energy is less 
than that required for the FS to explain the rest of the 
observations, i.e. the incoming ejecta provide only the fresh 
electrons to radiate at radio frequencies, but do not alter 
the dynamics of the FS. 

After including the 1 day radio flare measurements, as 
well as the millimetre measurements at 1-10 days, we also 
find with the Jet model a higher density solution (Fig. 4), 
n 1cm -3 , which is close to that obtained for other af- 
terglows. This higher density^is required by the larger FS 
peak flux necessary to accommodate the 1 day flare, when 
Fsghz = 0.36 mjy. For a uniform CBM, the high density 
fit is slightly better than the low density solution, the im- 
provement being more substantial for a wind. Note, however, 
that the former underestimates the radio measurements at 
30 days. 

The SO model yields fits of the same quality as the Jet 
model (Table 1). The best fit structural parameter satisfies 
q > q(p) for a uniform medium and q < q(p) for a wind (Ta- 
ble 3), thus the post-break afterglow emission arises mostly 



» It also lowers the cooling frequency, u c , to just slightly above 
the (blueward of) optical domain, which requires a hard electron 
distribution with index p it 2/3 ox ~ 4/3. Furthermore, a low 
electron parameter £; is required for the FS peak frequency to 
reach 10 GHz as early as 1 day after the burst, and to explain the 
radio flare 



from the outflow core in the former case and from the enve- 
lope for the latter. 

The best fit obtained with the EI model for a wind 
medium is shown in Fig. 5. x 2_w i se > this l° w density model 
provides a better fit than the Jet and SO models, although 
it overestimates the early radio measurements. Also, this so- 
lution has an extremely low wind density (Table 4), about 
100 times smaller than known for Galactic WR stars. For 
such a tenuous wind, equation (5) implies that, if the out- 
flow opening 6j et is wider than 2°, then the jet-break would 
appear later than about 30 days (i.e. after the latest mea- 
surement), and the corresponding outflow energy would be 
larger than 2 x 10 50 ergs, which is less than that of the jet 
shown in Fig. 4. The best fit obtained for a uniform medium 
is very poor (it either overestimates the radio emission or 
underproduces X-rays) and is not shown. 

3.3 990510 

The best with the Jet model for a uniform CBM is shown 
in Panaitescu & Kumar (2001). A wind provides a poorer 
fit (Table 1), as it fails to accommodate the strong break 
of the optical light-curves of this afterglow. Structured out- 
flows yield poorer fits for either type of medium, the best fit 
structural parameter being q = 1.8 for a uniform CBM and 
q = 1.4 for a wind. 

3.4 991216 

The best fit obtained with the Jet model and a uniform 
CBM is shown in Panaitescu & Kumar (2001). Because the 
radio decay of this afterglow, Fsghz oc £-°- 77 ± 06 ; j s s jg_ 
nificantly slower than that measured at optical wavelengths 
after the break, F a oc £ _1 - 66 ± 012 ; i n our previous modelling 
we have employed a double power-law electron distribution, 
harder (smaller index p) at low energies, to accommodate 
the shallower radio decay, and a softer (larger index p) at 
high energies, to explain the post-break optical light-curve 
decay. 

Here we use a single power-law electron distribution, to 
assess the ability of each model to explain all the data with- 
out recourse to a light-curve decay steepening originating 
in a break in the electron distribution. Consequently, the 
Jet model (as well as any other single emission component 
model) will not accommodate the decays of both the radio 
and optical emissions of the afterglow 991216, as illustrated 
by the higher Xv given in Table 1 and the best fit shown in 
Fig. 6. 

The addition of emission from the RS which energizes 
the ejecta catching up with the FS at about 1 day improves 
the radio fit, as shown in Fig. 7. A slower decay radio light- 
curve is obtained as the FS emission overtakes that from the 
RS at about 10 days. As for the afterglow 990123 (Fig. 3), 
the energy injected is small enough that it does not alter 
the dynamics of the FS. However, the afterglow parameters 
of the Jet+EI model for 991216 are different than those ob- 
tained with the Jet model (Fig. 6) because we require now 
that the FS peaks in the radio at a later time and at a lower 
flux (~ 0.1 mJy)H. 

II This requires a tenuous CBM, which increases the cooling fre- 
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The SO model yields a best fit which is very similar to 
that shown in Fig. 6 for the Jet model. The smaller intensity- 
averaged source size resulting in the SO model leads to a 
larger amplitude for the interstellar scintillation and, thus, 
to a better x 2 -wise fit (Table 1), but this model also fails to 
accommodate the slow radio decay. The best fit structural 
parameter satisfies q q(p) (Table 3), i.e. the outflow en- 
velope emission dominates over that from the core after the 
optical break. 

3.5 000301c 

The best fit for the Jet model is discussed in Panaitescu 
(2001). As for the afterglow 991216, the slower radio decay 
of 000301c, F$ GHz oc t- 1A±0 - [i , and its faster optical fall-off, 
F a oc i- 2 - 8 ' 5 ± - 12 j prompted us to consider a double power- 
law electron distribution. Furthermore, the large break mag- 
nitude Aa = 1.8 ± 0.2 exhibited by the optical light-curve 
000301c (the largest in the entire sample - Table 1) exceeds 
that allowed by the Jet model for the electron distribution 
index p required by the pre-break decay index and optical 
SED slope (Panaitescu 2005). Such a strong break is also 
better explained if there is a contribution from the passage 
of a spectral break. However, since we want to test the three 
models for light-curve break without any contribution from 
another mechanism, we list in Table 1 the best fit obtained 
previously with a soft electron distribution. 

Fig. 8 shows the best fit obtained with the SO model, 
which is only slightly better than that of the Jet model, has 
similar parameters, and shares the same deficiencies: it fails 
to explain the sharpness of the optical light-curve break and 
yields a radio emission decaying faster than observed. We 
note that optical measurements between 3.0 and 4.5 days, 
when a light-curve bump is seen, have been excluded from 
the fit. 

3.6 000926 

The best fit obtained with the Jet model is shown in 
Panaitescu & Kumar (2002). After including in the data set 
the near-infrared measurements, the noisy /-band, and some 
optical outliers which were previously excluded, the best fit 
is x 2 -wise poorer but the fit parameters are the same as be- 
fore (within their uncertainty). An equally good fit can be 
obtained with the SO model (Fig. 9). We note that for both 
the Jet and SO models, the X-ray emission is mostly inverse 
Compton scatterings. 

3.7 010222 

The best fit obtained with the Jet model is shown in Fig. 
10. Note that the model radio light-curve decays after 1 day 

quency (y c oc n _1 before the jet-break and v c oc n " 5 / 6 after 
that). Going from the Jet fit to the Jet+EI fit increases the cool- 
ing frequency from about 10 15 Hz to 3 X 10 17 Hz. Then, to ac- 
commodate the optical-to-X-ray SED slope, the Jet+EI model 
requires a softer electron distribution than the Jet model, which 
in turn yields a faster post-break decay of the afterglow light- 
curve, as shown in Fig. 7, and a poorer fit to the optical emission 
after 10 days 



faster than observed. For a uniform CBM, the best fit ob- 
tained with a structured outflow is slightly better in the 
radio, but it too fails to explain the slow radio decay. For 
a wind, the best fit obtained with the SO model overesti- 
mates the millimeter emission [F25oghz(1 — 100 d) ~ 1.2 
mJy, .F350 GHz (0.3 — 20 <i) ~ 3.9 mJy] attributed to the host 
galaxy by Frail et al. (2002), who also suggest a possible 
host synchrotron emission at 8.5 GHz of ~ 20 /xJy. This 
is only marginally consistent with the -Fi. 4 GHz (447 d) = 
— 1 ± 35/iJy reported by Galama et al. (2003) and the 
host synchrotron SED, F v oc v~ ' 75 , adopted by Frail et 
al. (2002), which imply a host flux Fsghz = ± 9 piJy. 

If we subtract a radio host contribution of Fv hoat ^ — 
20 {v/8 GHz)~ ' 75 /ijy from the 44 radio measurements, 
then, within the Jet model, the contribution to \ 2 of the 
radio data decreases from 107 to 72 for a uniform medium, 
and from 103 to 53 for a wind, with similar changes for 
the SO model. Though these are significant improvements, 
the best fit still underestimates the radio flux measured 
after 10 days, because the optical measurements have a 
smaller uncertainty and determine the electron index p and, 
implicitly, the decay of the model radio emission. In gen- 
eral, one-component models cannot accommodate both the 
host-subtracted radio emission, Fsghz oc ^-°- 76 ±°- 12 ] Q f 
the afterglow 010222 at 1-200 days, and the optical decay, 
F oc t' 1 - 78 * 08 , observed at 10-100 days. 

A better fit to the radio emission of the afterglow 010222 
may be obtained with the Jet model if, in addition to the 
FS emission, there is a RS contribution to the early radio 
afterglow, as illustrated in Fig. 11 for a uniform CBM. Just 
as for the afterglows 990123 and 991216, the incoming ejecta 
carry less energy than that in the FS (and do not alter its 
dynamics) and the medium density is lower than for the 
Jet model because the FS peak flux must be smaller, to 
match the radio flux measured after 10 days, when the FS 
synchrotron peak frequency crosses the radio domain. 



3.8 011211 

As illustrated in Fig. 12, even for a uniform CBM, the Jet 
model has difficulty in accommodating the sharp break ex- 
hibited by the optical light-curve of this afterglow, assum- 
ing that the 1-2 day optical emission is not a fluctuation. 
Such fluctuations have been seen in the afterglows 000301c, 
021004, and 030329, but the lack of a continuous monitoring 
for 011211 prevents us to determine if its optical emission at 
1 day is indeed a fluctuation. A stronger break can be ob- 
tained with the SO model when the brighter, outflow core 
becomes visible to an observer located outside it, hence the 
better fit that the SO model yields (Fig. 13). 



3.9 020813 

The best fit obtained with the Jet model is shown in Fig. 
14. Structured outflows yield fits with q > q(p) for a uni- 
form CBM and q < q(p) for a wind, i.e. with the post-break 
emission arising mostly in the outflow core and envelope, 
respectively. A significant improvement over the Jet model 
is obtained only for a uniform CBM (Table 1). 

A fit of comparable quality is also obtained with the 
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EI model for a wind medium (Fig. 15), but not for a uni- 
form CBM, which overproduces radio emission and largely 
underestimates the observed X-ray fluxes. For the best fit 
parameters given in Table 4, equation (5) shows that a jet- 
break would occur after the last measurement (40 days) if 
the outflow opening were larger than 5°, which implies a 
minimum collimated energy of 2 x 10° ergs, which is 3 times 
larger than that of the jet shown in Fig. 14. Therefore, just 
as for 990123, the EI model may not involve significantly 
larger energies than the Jet model, if allowance is made for 
a possible outflow collimation. 

3.10 030226 

As for 011211, the strong optical light-curve break of this 
afterglow is better accommodated by the SO model than 
the Jet model (Fig. 16 and 17), though an acceptable fit to 
the optical break is obtained only for a uniform CBM (in 
addition, for a wind medium, the SO model overproduces 
radio emission). However, given the poor coverage before 
0.3 days, it is possible that the measured optical emission 
at that time is a fluctuation below that expected in the Jet 
model. 



3.11 970508 

Although the decay of the optical emission of this afterglow 
does not exhibit a steepening, being so far the longest mon- 
itored power-law fall-off (2-300 days), it can be explained 
with the Jet model if its unusual rise at 1-2 days is inter- 
preted as a jet becoming visible to an observer located out- 
side the initial jet opening (see fig. 1 in Panaitescu & Kumar 
2002). This requires a rather wide jet, with 9j et = 18°, and 
yields a poor fit to the X-ray emission (xt = 2.8). A less 
energetic envelope surrounding this jet would produce the 
prompt 7-ray emission and would account for the flat opti- 
cal emission prior to the 2 day optical rise, hence this model 
is similar to a Structured Outflow with a wide uniform core 
and an envelope whose structure is poorly constrained by 
observations. 

An other possible reason for the 2 day sharp rise of this 
afterglow is a powerful energy injection episode in which the 
incoming ejecta carry an energy larger than that already 
existing in the FS (Panaitescu, Meszaros & Rees 1998), i.e. 
the blast-wave dynamics is altered drastically by the energy 
injection. Fig. 18 shows the best with the EI model. Note 
that the given that the light-curve rise is achromatic, as it 
arises from the blast-wave dynamics, which explains why the 
X-ray light-curve cannot be accommodated. 



4 CONCLUSIONS 

In this work we expand our previous investigation 
(Panaitescu & Kumar 2003) of the Jet model to four other 
GRB afterglows for which an optical light-curve break was 
observed: 010222 (for which radio observations were not 
available at the time of our previous modelling), 011211, 
020813, and 030226. We also note that, here, we have not 
employed a broken power-law electron distribution to ex- 
plain both the slower decaying radio emission and the faster 



falling-off optical light-curve observed for some afterglows 
(e.g. 991216 and 000301c). 

Our prior conclusion that a homogeneous circumburst 
medium provides a better fit to the broadband afterglow 
emission than environments with a wind-like stratification 
stands (Table 1). For afterglows with a sharp steepening Aa 
of the optical light-curve decay (990510, 000301c, 011211, 
and 030226), this is due to that the transition between the 
jet asymptotic dynamical regimes (spherical expansion at 
early times and lateral spreading afterward) occurs faster 
for a homogeneous medium than for a wind (Kumar & 
Panaitescu 2000). For the afterglows with shallower breaks, 
a wind medium (as expected if the GRB progenitor is a mas- 
sive star) provides an equally good fit as (or, at least, not 
much worse than) a homogeneous medium. 

For the acceptable fits (Table 2) obtained with the Jet 
model and for a uniform medium, the jet initial kinetic en- 
ergy, Ej e t, is between 2 and 6 xlO 50 ergs, the initial jet 
opening angle, 6j et , is between 2° and 3° (with 000926 an 
outlier at 8°), and the medium density, n, between 0.05 and 
1 cm~ 3 (with 000926 an outlier, again, at 20 cm -3 ). For 
a wind, Ej et has about the same range, 9j et ranges from 
2° to 6°, while the wind parameter A* is between 0.1 and 
2. 70% of the 64 Galactic WR stars analyzed by Nugis & 
Lamers (2000) have a A* parameter in this range; for the 
rest A* £ (2,5). We note that the Ej e t resulting from our 
numerical fits is uncertain by a factor 2, 9j e t by about 30%, 
n by almost one order of magnitude, and A„ by a factor of 
2. 

Although the inferred wind density parameter A* is 
in accord with the measurements for WR stars a uniform 
medium is instead favoured by the fits to the afterglow data. 
This is an interesting issue for the established origin of long 
bursts in the collapse of massive stars. Wijers (2001) has 
proposed that the termination shock resulting from the in- 
teraction between the wind and the circumstellar medium 
could homogenize the wind. The numerical hydrodynamical 
calculations of Garcia-Segura, Langer & Mac Low (1996) 
show that, for a negligible pressure nT of the circumstel- 
lar medium (about 10 4 cm _3 K), the termination shock of 
the RSG wind has a radius of 10 pc, outside a shell of uni- 
form density is form. However, this radius is much larger 
than the distance of about 1 pc where the afterglow emis- 
sion is produced. Chevalier, Li & Fransson (2004) made the 
case that, if the burst occurred in an intense starburst re- 
gion, where the interstellar pressure is about 10 7 cm _3 K, 
then the shocked RSG wind would form bubble of uni- 
form density of about 1 cm~ 3 extending from 0.4 to 1.6 
pc, in accord with the fit densities and the blast-wave ra- 
dius r = 8cr 2 t/(l + z) = 0.5 (E , 53 /no) 1/4 [t d /(l + z)] 1/4 pc, 
resulting from equation (2). 

The Structured Outflow model involving an ener- 
getic, uniform core, surrounded by a power-law envelope 
(equation [7]) which impedes the lateral spreading of the 
core, retains most of the ability of the Jet model to yield 
a steeply decaying optical light-curve after the core (or its 
boundary) becomes visible to an observer located outside 
the core opening (or within it), while accommodating, at 
the same time, the observed slope of the optical spectral 
energy distribution (SED). This is so because, in both mod- 
els, more than half of the increase Aa in the light-curve 
power-law decay index is due to the finite opening of the 
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bright(er) ejecta. For a jet, the remainder of the steepening 
Aa is caused by the lateral spreading. For a structured out- 
flow, an extra contribution to the steepening, which can lead 
to a Aa even larger than that produced by a jet, results if 
the observer is located outside the core opening. 

For a uniform medium, we find that the Structured Out- 
flow model provides a better fit than the Jet model for 6 
of the 10 afterglows analyzed here (Table 1). For a wind, 
the former model works better for three afterglows, fits of 
comparable quality being obtained for the other six. Only 
the afterglow 990510 is better explained with a jet, for ei- 
ther type of circumburst medium. Within the framework of 
structured outflows, 8 out of 10 afterglows are better fit with 
a uniform medium than with a wind, which is partly due to 
that, just as for a jet, the light-curve steepening is sharper 
if the medium is uniform (Panaitescu & Kumar 2003). 

For the "acceptable" fits (xt < 4) obtained with the 
SO model, the ejecta kinetic energy per solid angle in the 
core, Eo, ranges from 0.5 to 3 x 10 53 ergs sr _1 for either type 
of medium, with an uncertainty of a factor 3. The angular 
opening of the core, 6 core , is between 0.5° and 1° for a ho- 
mogeneous medium, and about 1° for a wind (991216 be- 
ing an outlier at 0.5°), with an uncertainty less than 50%. 
The structural parameter q (equation [7]), which character- 
izes the energy distribution in the outflow envelope, is found 
to be between 1.5 and 2.7 for a uniform medium, and be- 
tween 1.0 and 2.3 for a wind (with 000926 being an outlier 
at 2.9). These values are consistent with that proposed by 
Rossi et al. (2002), q = 2, to explain the quasi-universal 
jet energy resulting in the Jet model. For structured out- 
flows, the ejecta kinetic energy contained within an opening 
of twice the location of the observer (relative to the outflow 
symmetry axis), 9 obs , ranges from 1 to 4 xlO ergs, for ei- 
ther type of medium. Hence, the energy budget required by 
the Structured Outflow model is very similar to that of the 
Jet model. The best fit medium densities obtained with the 
Structured Outflow model are similar to those resulting for 
a Jet. 

We find that the Energy Injection model, where a 
light-curve decay steepening is attributed to the cessation of 
energy injection in the forward shock, can be at work only in 
two afterglows, 990123 and 020813. The major shortcoming 
of the EI model is that, in order to explain the steep fall-off 
of the optical light-curve after the break with a spherical 
outflow, an electron distribution that is too soft (i.e. an in- 
dex p too large) is required to accommodate the relative 
intensities of the radio and X-ray emissions. Low density 
solutions [n < 1 cm~ 3 , A„ < 0.1), for which the X-ray emis- 
sion is mostly synchrotron, yield a cooling frequency that 
is too low, resulting in model X-ray fluxes underestimating 
the observations by a factor of at least 10. High density solu- 
tions (n > 10cm~ 3 , A* > 1), for which the X-ray is mostly 
inverse Compton scatterings, produce a forward shock peak 
flux that is too large, which leads to model radio fluxes over- 
estimating the data by a factor of 10. 

We note that, due to adiabatic cooling, the GRB ejecta 
electrons which were energized by the reverse shock during 
the burst, do not contribute significantly to the afterglow 
radio emission at days after the burst. Numerically, we find 
that, if the microphysical parameters are the same for both 
the reverse and forward shocks, then the 1 day radio flare 
of the afterglow 990123 cannot be explained with the emis- 



sion from the GRB ejecta electrons, as proposed by Sari & 
Piran (1999), as the synchrotron emission from these elec- 
trons peaks in the radio earlier than 0.1 day, regardless of 
how relativistic is the reverse shock. For the reverse shock 
emission to be significant at 1 day after the burst, there must 
be some other ejecta catching up with the decelerating out- 
flow at that time or the assumption of equal microphysical 
parameters for both shocks must be invalid (Panaitescu & 
Kumar 2004). Furthermore, the optical emission from the 
ejecta electrons is below that observed at 100 s (Akerlof et 
al. 1999) by a factor 5 for a wind and a factor 50 for a 
uniform medium. 

A delayed injection of ejecta, carrying less kinetic en- 
ergy than that already existing in the forward shock (i.e. 
not altering the forward shock dynamics), may be at work 
in the afterglows 990123, 991216, and 010222, whose radio 
decay is substantially slower than that observed in the op- 
tical after the break. As shown by us (Panaitescu & Kumar 
2004), such a decoupling of the radio and optical emissions 
cannot be explained by models based on a single emission 
component (forward shock) and require the contribution of 
the reverse shock. Figs. 3, 7, and 11, illustrate that the sum 
of the reverse and forward shock emissions yields a slower 
decaying radio emission, consistent with the observations. 
However, for the forward shock peak flux to be as dim as 
observed in the radio after 10 days, when the forward shock 
peak frequency decreases to 10 GHz, an extremely tenuous 
medium is required: n ~ 10~ 3 cm~ 3 or A, ~ 0.05. The for- 
mer might be compatible with a massive star progenitor if 
GRBs occur in the "superbubble" blown by preceding su- 
pernovae exploding in the same molecular cloud (Scalo & 
Wheeler 2001), while the latter points to a progenitor with 
a low mass and low metallicity or, otherwise, to a small mass- 
loss rate and high wind speed during the last few thousand 
years before the collapse, when the environment within 1 pc 
is shaped by the stellar wind. 

As a final conclusion, we note that the Structured Out- 
flow model is a serious contender to the Jet model in accom- 
modating the broadband emission of GRB afterglows with 
optical light-curve breaks. In both models, the best fit pa- 
rameters describing the ejecta kinetic energy - jet opening 
& energy or core opening & energy density along the out- 
flow symmetry axis - have narrow distributions, hinting to 
a possible universality of these parameters. 
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Figure 1. Best fit for the afterglow 980519 with the SO model. 
The data set contains measurements at 1.4, 4.9, 8.5, 100 GHz 
(Frail ct al. 2000b), /, R, V, B, U bands (Halpern et al. 1999, Vrba 
et al. 2000, Jaunscn ct al. 2001), and 5 keV (inferred from the 2- 
10 keV fluxes presented by Nicastro et al. 1999). Vertical bars 
on the model radio light-curve indicate the amplitude of Galactic 
interstellar scintillation. The model X-ray emission is mostly in- 
verse Compton scatterings. For this afterglow, a redshift was not 
measured. We assumed z = 1. 



Figure 2. Best fit for the afterglow 990123 with the Jet model 
and emission from the ejecta electrons accelerated by the reverse 
shock during the burst (t ~ 100 s). The data set contains mea- 
surements at 1.4, 4.9, 8.3, 15, 86, 353 GHz (Galama et al. 1999, 
Kulkarni et al. 1999), K, H, I, R, V, B (Akcrlof ct al. 1999, 
Galama et al. 1999, Kulkarni ct al. 1999), and 6 keV (inferred 
from the 2-10 keV fluxes reported by Piro 2000). Triangles de- 
note 2(7 observational limits on the 8.3 GHz emission. Radio fluxes 
have been shifted upwards by a factor 100, for clarity. For this 
model, the ejecta Lorentz factor Yi was not constrained by equa- 
tion (13), because the assumption of instantaneous ejecta release 
may not be valid, and was left a free parameter. The optical 
peak at 100 s and the radio flux before 1 day represent the re- 
verse shock emission, which is later overshincd by the forward 



shock. Model parameters are: Eo = 3 X 10 : 
n = 0.004cm- 3 , E, = 1.3 X 10 53 crgs sr" 1 



50 ergs sr 1 , 



7 jet 



2.2° 



r\ = 5,000, t off = 90 



s, e B = 4 X 10~ 4 , Si = 0.075, p = 2.2, A v = 0.08 for a uniform 
medium and E = 1.5 X 10 50 crgs sr -1 , 6 jet = 6.0°, A„ = 0.06, 
F t = 3 x 10 53 crgs sr" 1 , T; = 1,200, t off = 170 s, e B = 7x 10" 5 , 
Ei = 0.080, p = 2.1, A v = 0.07 for a wind. The addition of the 
early optical data has worsened the fits: xi = 14 for the former 
medium, xl = H f° r tne latter. Note that the wind model ac- 
commodates most of the early optical observations but, for either 
type of medium, the reverse shock emission peaks in the radio too 
early, at about 0.1 day. 
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Figure 3. Best fit for the afterglow 990123 with the Jet+EI 
model for a uniform medium and a wind. The peak of the radio 
light-curve before 1 day represents the reverse shock emission, 
while the forward shock peak frequency crosses the radio domain 
after 10 days. 




time (day) 

Figure 4. Best fit for the afterglow 990123 with the Jet model. 
Note that, in contrast to the fits shown in figure 3, the medium 
density of these fits is similar to that determined for other af- 
terglows. However, the electron energy distribution has unusual 
parameters and p. 



Figure 5. Best fit for the afterglow 990123 with the EI model for 
a wind medium. The model radio emission until 10 days arises in 
the reverse shock which has energized the incoming ejecta, after 
which the forward shock emission overtakes it. 
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Figure 6. Best fit for the afterglow 991216 with the Jet model. 
The data set consists of measurements at 1.4, 4.9, 8.5, 15 GHz 
(Frail ct al. 2000c), K, H, J, I, R, V bands (Garnavich ct al. 2000, 
Halpcrn ct al. 2000), and 6 keV (inferred from the 2-10 keV fluxes 
reported by Corbet & Smith 1999, Piro ct al. 1999, Takcshima ct 
al. 1999). Note that this model cannot explain the slower radio 
decay observed at about the same time with the faster optical 
fall-off. 
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Figure 7. Best fit for the afterglow 991216 with the Jct+EI 
model for a uniform medium and a wind. The addition of the 
reverse shock, which yields the 2 day peak of the radio emission, 
improves the fit to the radio data compared to that obtained with 
the Jet model (figure 6). 
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Figure 9. Best fit to the afterglow 000926 with the SO model. 
The data consists of measurements at 1.4, 4.9, 8.5, 15, 22, 99 
GHz (Harrison et al. 2001, Frail ct al. 2003), K, H, J, I, R, V, B, U 
bands (Fynbo ct al. 2001, Harrison ct al. 2001, Price et al. 2001), 
0.5 and 3 keV (inferred from the 0.2-1.5 keV and 2-10 kcV fluxes 
reported by Harrison et al. 2001 and Piro ct al. 2001). 



10' 
10° 
10" 1 

Si 

~* 2 

JL 10 

> 

LL 

10" 3 
10" 4 
10" F 



8.5 GHz 



1 1 1 
unif CBM= 

wind 




0301. jo 



10° 



10 1 10 2 

time (day) 



Figure 8. Best fit for the afterglow 000301c for the SO model. 
The data set contains measurements at 1.4, 4.9, 8.5, 15, 22, 100, 
250 350 GHz (Berger ct al. 2000, Smith ct al. 2001, Frail ct al. 
2003), and K, J, I, R, V, B, U bands (Jensen et al. 2001, Rhoads 
& Fruchter 2001). The 250 GHz fluxes and their 2cr limits (tri- 
angles) have been shifted downward by a factor 20, for clarity. 
Measurements between 3.0 and 4.5 days, when there is a fluc- 
tuation in the optical emission, have been left out. The best fit 
parameters for the uniform medium are given in Table 3, those 
for the wind medium are Eo = 1.5 X 10 53 ergs sr~ 1 , Score = 0.9°, 
dobs/dcore = 2.6, q = 1.4, A, = 0.6, £ B =3x 10~ 3 , e t = 0.018, 
p = 2.4, A v = 0.06. 
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Figure 10. Best fit to the afterglow 010222 with the Jet model. 
The data set contains measurements at 1.4, 4.9, 8.5, 15, 22 GHz 
(Galama et al. 2003), K, J, I, R, V, B, U bands (Masetti et al. 
2001, Stanek et al. 2001), and 5 kcV (inferred from the 2-10 keV 
fluxes reported by in't Zand et al. 2001). The steep X-ray fall-off 
after 2 days for a wind medium is due to a high energy cutoff of 
the electron distribution, corresponding to a total electron energy 
of 50 per cent of the post-shock energy. 
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Figure 11. Best fit for the afterglow 010222 with the Jet+EI 
model. For a uniform medium the forward shock component ac- 
commodates better the radio data after 10 days than for a wind. 
The 1 day hump in the radio light-curve is the reverse shock emis- 
sion, the forward shock overtaking at 10 days. 



Figure 13. Best fit for the afterglow 011211 with the SO model. 
The best fit parameters for the uniform medium are given in 
Table 3, those for the wind medium are Eq = 1.2 X 10 53 crgs sr -1 , 
C ore = 0.9°, d obs /e core = 3.0, q = 2.3, A* = 0.7, e B = lO" 3 , 
Si = 0.0081, p = 2.3, A v = 0. 




Figure 12. Best fit for the afterglow 011211 with the Jet 
model. The data set contains measurements at 8.5, 22 GHz 
(http://www.aoc.nrao.cdu~dfraiV011211.dat), K, J, I, R, V, B, U 
bands (Holland et al. 2002, Jakobsson et al. 2003), and 2 keV 
(inferred from the 0.2-5 keV fluxes measured by Borozdin & 
Trodolyubov 2003). The best fit parameters are Eo = 2 X 
10 52 crgs sr" 1 , 9 jet = 5.5°, n = LOcm" 3 , e B = 6 X 10" 4 , 
£i = 0.047, p = 2.3, Ay = for a uniform medium and 
E = 4 x 10 52 ergs sr" 1 , 9 jet = 4.0°, A, = 0.6, e B = 5 x 10~ 4 , 
£i = 0.030, p = 2.2, Ay =0 for a wind. 



Figure 14. Best fit for the afterglow 020813 with the 
Jet model. The data set contains measurements at 4.5, 
8.5 (http://www.aoc.nrao.cdu~dfrail/grb020813.dat), 250 GHz 
(Bertoldi et al. 2002), K, H, J, I, R,V, B,U bands (Covino et al. 
2003, Gorosabel et al. 2004), and 2 keV (inferred from the 0.6-6 
keV fluxes measured by Butler et al. 2003). For the latter, the 
sharp steepening of the X-ray light-curve at <; 1 day is due to 
the electron high-energy cut-off. 
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Figure 15. Best fit for the afterglow 020813 with the EI model 
for a wind medium. The reverse and forward shock radio emis- 
sions are equal at 1 day, after that the forward shock emission is 
stronger. 



Figure 17. Best fit for the afterglow 030226 with the SO model. 
The best fit parameters for the uniform medium are given in 
Table 3, those for the wind medium are Eq = 1.3 X 10 53 crgs sr -1 , 
9 core = 0.5°, e obs /9 core = 4.1, q = 2.8, A* = 0.5, e B = 6 x IO" 4 , 
Si = 0.017, p = 2.2, A v = 0.03. 
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Figure 16. Best fit for the afterglow 030226 with the 
Jet model. The data set contains measurements at 8.5, 
15, 22 GHz (http://www.aoc.nrao.cdu~dfraiVgrb030226.dat), 
K,H,J,I,R,V,B,U bands (Klose et al. 2004, Pandey ct al. 
2004), and 2 keV (inferred from the 0.3-10 keV count rate re- 
ported by Klose et al. 2004). The best fit parameters are Eq = 
6 X 10 52 crgs sr -1 , 6 jet = 2.8°, n = 1.6cm" 3 , e B = 2 X 10~ 4 , 
Ei = 0.029, p = 2.2, Ay = 0.09 for a uniform medium and 
E = 9 x 10 52 ergs sr" 1 , 9 jet = 2.6°, A, = 0.5, e B = 4 x 10~ 4 , 
Ei = 0.016, p = 2.1, A v = 0.04 for a wind. 
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Figure 18. Best fit for the afterglow 970508 with the EI model. 
The data set contains 309 measurements at 1.4, 4.9 GHz (Frail, 
Waxman & Kulkarni 2000a), 86 GHz (Bremer et al. 1998), 
K,I, R, V, B bands (Chary et al. 1998, Galama ct al. 1998, Gar- 
cia ct al. 1998, Pedersen ct al. 1998, Sahu et al. 1997, Sokolov ct 
al. 1998), and 5 keV (inferred from 2—10 keV fluxes and spec- 
tral slope given in Piro et al. 1998). The fit parameters arc: 
E % = 9 X 10 51 ergs sr^ 1 , e = 2.4, t off = 1.6 d, n = 0.9cm- 3 , 
£b = 

ixl = 4.5) and E t 
A, = 0.7cm- 3 , E B = 2 x 10 
for a wind (x 2 , = 5.9). 
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